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I. INTRODUCTION 


Just after ~ 10“^^ seconds of the big bang, the universe was hlled with a state of matter 
known as quark-gluon plasma. The quark-gluon plasma is much hotter and denser than 
the ordinary matter we see today. The temperature of the quark-gluon plasma is ^ 200 
MeV 3.2 X 10^^ K) which is about million times larger than the temperature of the 
sun (~ 1.56 X 10^ K). The quark-gluon plasma is much denser than neutron stars, besides 
black holes, there’s nothing denser than this. Hence it is important to recreate this early 
universe scenario in the laboratory, i. e., to produce quark-gluon plasma in the laboratory. 
At present RHIC (relativistic heavy-ion colliders) at BNL and LHC (large hadron collider) 
at CERN are the best facilities to produce quark-gluon plasma in the laboratory 1, 

RHIC and LHC use high energy heavy-ion collisions to produce quark-gluon plasma. 
RHIC collides two gold nuclei at ^/s = 200 GeV per nucleon j^, 4] with total energy ~ 40 
TeV. The LHC (in its hrst run) collides two lead nuclei at ^/s = 2.76 TeV per nucleon M 
with total energy ~ 574 TeV. Since these huge total energies are deposited in very small 
volume at the initial moment of the nuclear collisions at RHIC and LHC, the energy density 
necessary to create temperature ~ 200 MeV to produce quark-gluon plasma might have 
been achieved at RHIC and LHC. The LHC (in its second run) will collide two lead nuclei 
at ^/s = 5.5 TeV per nucleon with total energy ~ 1150 TeV which will be even higher. 

However, since the two nuclei at RHIC and LHC travel almost at speed of light, the quark- 
gluon plasma produced at RHIC and LHC may be in non-equilibrium. This is because the 
partons inside the nuclei at RHIC and LHC carry very high longitudinal momentum and very 
small transverse momentum just before the nuclear collision. Since the hadronization takes 
place in a very short time scale (the typical time scale for hadronization in QCD is ~ 10“^^ 
seconds), there may not be enough time for these highly non-isotropic partons to thermalize 
at RHIC and LHC. Hence, in order to make meaningful comparison of the theory with 
the experimental data on hadron production it may be necessary to study nonequilibrium- 
nonperturbative QCD at RHIC and LHC. This, however, is a difficult problem. 

Since quarks and gluons are not directly observed we can not directly detect the quark- 
gluon plasma at RHIC and LHC. Hence indirect signatures are proposed for the detection of 
quark-gluon plasma at RHIC and LHC. The main signatures of quark-gluon plasma detection 
are 1) j/V’ suppression/production 2) dilepton production 3) direct photon production, 4) 
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strangeness enhancement and 5) jet quenching etc.. In this paper we will focus on the j/'ip 
suppression/production signature for the quark-gluon plasma detection at RHIC and LHC. 

Note that the i/pj suppression was predicted to be a signature of quark-gluon plasma 
detection by using lattice QCD calculation at finite temperature in equilibrium j^. It was 
argued that at high temperature the Debye screening length is much smaller than the j/pj 
radius, leading to complete suppression oi j/p) production in quark-gluon plasma. However, 
this prediction uses lattice QCD calculation at finite temperature which may not be appli¬ 
cable at RHIC and LHC where the quark-gluon plasma may be in non-equilibrium. Hence, 
in order to make comparison of the theory with the experimental data on j/pj production 
at RHIC and LHC, it may be necessary to study j/V’ production in non-equilibrium QCD. 

The suppression is experimentally observed in An-An collisions at i/i = 200 GeV at 


RHIC jo,!!^ and in Pb-Pb collisions at ^/s = 2.76 TeV at LHC (in its first run) llll-ll3|. It is 


interesting to note that in most central heavy-ion collisions the PHENIX collaboration [9[ at 
RHIC has measured more j/ip suppression than the ALICE collaboration at LHC (in its 
first run). In order to explain this it is argued that there is enhancement [l^ and charm 
recombination in QCD medium etc., however these approaches need to incorporate non- 
perturbative QCD mechanism for heavy quarkonium production from cc pair. Hence the 
determination of in-medium properties of heavy quarkonium at RHIC and LHC remains a 
challenging theoretical task. 

The experimental data at RHIC and LHC suggest that the j/ip production cross section 
is modified in the heavy-ion collisions in comparison to the p-p collisions at the same center 
of mass energy. This modification can be due to the cold nuclear matter effects such as 
nuclear shadowing, saturation, energy loss and quark-antiquark breakup etc. and due to the 
hot QCD medium effects such as color screening and secondary charm production etc.. The 
cold nuclear matter effects are studied both theoretically and experimentally, for example 
by using p-A collisions at the corresponding center of mass energy. The hot QCD medium 
effects such as color screening is well studied theoretically by using lattice QCD at finite 
temperature js, 17|, but needs to be extended to non-equilibrium QCD to be applicable at 
RHIC and LHC. Similarly secondary charm production is well studied at finite temperature 
14j. IlSl. IlSl. [lOj, but how these secondary charm-anticharm produces j/ip by using non- 
perturbative QCD mechanism is not well understood, which also needs to be extended to 
non-equilibrium QCD at RHIC and LHC. 
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Hence the j/'ip production mechanism in QCD in vacuum (for example in p-p collisions) 
need to be modihed to include QCD medium effects at the heavy-ion collisions at RHIC and 
LHC. As mentioned above, since the two nuclei at RHIC and LHC travel almost at speed 
of light, the quark-gluon plasma at RHIC and LHC may be in non-equilibrium. Hence, in 
order to detect quark-gluon plasma by using as a signature, one needs to study i/pj 
production in non-equilibrium QCD at RHIC and LHC. 

Note that in order to study j/^p production in p-p collisions at high energy one needs to 
prove factorization theorem. The need to prove factorization theorem arises because in the 
absence of the proof of factorization theorem one will predict inhnite cross section of j/pj 
production 2CH28| . This is due to the interaction between charm (anticharm) quark with 
the nearby light quark or gluon. The soft gluon exchange between charm (anticharm) quark 
and the nearby light quark (or gluon) gives infrared divergence which makes the partonic 
level cross section inhnite. Hence it is important to prove that the effect of such infrared 
(soft) gluon exchanges either cancel or the non-canceling infrared divergences be absorbed 
into the dehnition of the non-perturbative matrix element of the j/p) because (soft) infrared 
regime corresponds to non-perturbative QCD. 

In the non-equilibrium QCD medium at RHIC and LHC there are many more nearby 
light quarks and gluons (in comparison to the corresponding situation in p-p collisions in 
QCD in vacuum) with which the charm (anticharm) quark interacts. Because of this, in¬ 
frared divergences seem to be severe at RHIC and LHC heavy-ion collisions than that in 
corresponding p-p collisions. Hence it is necessary to prove that the infrared divergences 
due to the soft gluons exchange between charm (anticharm) quark and the nearby light 
quark (or gluon) in non-equilibrium QCD medium at RHIC and LHC either cancel or the 
non-canceling infrared divergences be absorbed into the dehnition of the non-perturbative 
matrix element of the j/p) production in non-equilibrium QCD. 


Since the j/V’ suppression as a signa 
using j/p) from color singlet cc pair in 


iture of quark-gluon plasma detection was suggested by 
81], by comparing with the correlation length obtained 
from lattice QCD 29|, we will consider j/ip production from color singlet cc pair in non¬ 
equilibrium QCD in this paper. We will prove, in this paper, that the infrared divergences 
due to the soft gluons exchange between charm quark and the nearby light-like quark (or 
gluon) in non-equilibrium QCD at RHIC and LHC exactly cancel with the corresponding 
infrared divergences due to the soft gluons exchange between anticharm quark and the same 
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nearby light-like quark (or gluon) at all order in coupling constant in the j/'ip production 
from color singlet cc pair in non-equilibrium QCD. This proves the factorization theorem of 
j/p) production in non-equilibrium QCD at RHIC and LHC at all order in coupling constant. 

The paper is organized as follows. In section II we briefly discuss path integral formulation 
of non-equilibrium QCD. In section III we discuss the non-perturbative matrix element of 
j/p) production from color singlet charm-anticharm pair in non-equilibrium QCD. In section 
IV we discuss how longitudinal polarization of the gauge held in quantum held theory is 
used to describe infrared divergence. In section V we discuss the gauge held generated by 
the eikonal current of the light-like Wilson line in quantum held theory. In section VI we 
prove factorization theorem of j/p) production at high energy colliders. In section VII we 
prove factorization theorem of j/p) production in non-equilibrium QCD at RHIC and LHC. 
Section VHI contains conclusions. 

II. NON-EQUILIBRIUM QCD USING CLOSED-TIME PATH INTEGRAL FOR¬ 
MALISM 

Unlike pp collisions, the ground state at RHIC and LHC heavy-ion collisions (due to the 
presence of a QCD medium at initial time t = tin (say tjn=0) is not a vacuum state |0 > 
any more. We denote \in > as the initial state of the non-equilibrium QCD medium at tin- 

Consider massless scalar held theory hrst. In the closed-time path (CTP) formalism 
the generating functional in the path integral formulation for scalar held theory in non- 
equilibrium is given by 

Z[p, J_|_, J_] = J [d(j)^][d(j)_] exp[i[5'[0+] — >S'[0_] -1- J d^xJ+cp^ — J d^xJ_(f)_]] < 0+,O|p|O,0_ > 

where p is the initial density of state, S'[0] is the full action in scaler held theory and |0±, 0 > 
is the quantum state corresponding to the held conhguration (t>±{x,t = 0). Since there are 
two sources J_|_ and J_ corresponding to two time branches -|- and — in non-equilibrium 
quantum held theory one hnds that there are four Green’s functions in non-equilibrium. The 
four Green’s functions in non-equilibrium can be obtained from the generating functional as 
follows 

G++{x,x') = =< in\T(j){x)(P{x')\tn > 
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G _ (x, x') 

G+-{x,x') 

G-+{x,x') 


{-iy5j.{x)j.{x') 

5Z[p,J^,J_\ 

—i^5J+{x)J-{x') 

5Z[p,J^,J_] 

—i^5J_{x)J+{x') 


=< m|T0(x)0(a;')|m > 

< m|0(a;')0(a;)|m > 

< m|0(a;)0(x')|m > 


( 2 ) 


where +(— ) sign corresponds to upper(lower) time branch of the Schwinger-Keldysh closed¬ 


time path 
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J, l52(], T is the time order product and T is the anti-time order product. 

The generating functional in non-equilibrium QCD (including heavy quark) in the path 
integral formulation is given by 53|, [54 1 


Z[p, J+, J_, Pi, pt V-, vl, c, pt, 5+, r] 


= j lllQ+l MQ-] Ml/'l-l |#1+] 1#1-] l<ii’2+] ldi’2-] 1#2+] l<i<l>2-] ldi’3+] 1#3-] Ml/'S+I 

exp[ 3 1 d-il-icFtxlC+l - F“;;„[Q-]) - T(( 5 „Q 7)2 - 

3 3 

Z=1 Z=1 

+^4tYd^ -M + z?TVQ“+]^+ - -M + z7TVQ“-]^- + J+ ■ Q+ - J- ■ Q- 

3 

+ Y^[pi+ ■ '4)1+ - pi_ ■ 4)i_ + 4)1+ ■ pi+ - 4)i_ ■ pi_] -F ph+ ■ T+ - ph- • f + • ph+ - f _ ■ ph_]] 

z=i 


X < g+,V'l^1,V'l, ^+, ^+, oi p io, T_, T_, ^1,^1, r_,r-, q- > 


(3) 


where p is the initial density of state. The state |g^, "0^, 0 > corresponds to the field 

conhgurations g“(x,f = tin = 0), '4){x,t = tin = 0) and '4){x,t = tin = 0) respectively and 
J^°‘{x) is the external source for the quantum gluon field Q^°‘{x) and the p^, p'^, p^ are 
external sources for / = 1, 2, 3 = u, d, s quark fields respectively and ph is the external source 
for the heavy quark held T, mi is the mass of the light quark of type I, the M is the mass 
of the heavy quark and 


FLIQ] = 8 ,<rM - dM(x) + ar^QUxlQlix), 


pa2 


flU 


[01 = F»“'“[0]F“ [0], (4) 


Note that we work in the frozen ghost formalism 53|, l5^ for the medium part at the initial 


time t = tin = 0. 
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For the heavy quark Dirac held \l/r(a;) in non-equilibrium QCD, the nonequilibrium- 
nonperturbative matrix element of the type < m|^r(a;i)\l/r(a^i)'hs(a; 2 )tl/s(a; 2 )|m > is given 
by bl 


< m|\['r(xi)\l'r(xi)\['s(x2)tt's(x2)|m > 
6 6 6 6 


6r]hr{xi) 6r]hr{xi) 6r]hs{x2) 6r]hs{x2) 


Z[p, J+, J_, TJI, 77“ , C, pt pi, ptv% h-, vl, c 


5 V+^ V—1 V+1 V—\ I J+=J_=r7^=r;“=77“=r)“ =77^=r;|^=r7^=r;^=r;^=r;f_=r;^=T;f_=T7^=J7^=r;^j^=r;^^=0 


( 5 ) 


where the repeated (closed-time path) indices r, s = -|-, — are not summed and the suppres¬ 
sion of the normalization factor Z[0] is understood as it will cancel in the hnal result, see 
eq. 


III. NON-PERTURBATIVE MATRIX ELEMENT OF J/^ PRODUCTION FROM 
COLOR SINGLET CHARM-ANTICHARM PAIR IN NON-EQUILIBRIUM QCD 

The charm-anticharm quark (cc) pair production in the initial collision of the two nuclei 
at RHIC and LHC can be calculated by using pQCD. The spectral notation of j/t/' is given 
by where spin quantum number S = 1 (spin triplet), orbital angular momentum 

quantum number L = 0 (F-wave) and total angular momentum quantum number J = 1. 
Once the cc pair is formed in the color singlet and spin triplet state it can form j/'ip- The 
formation of j/p) from cc involves non-perturbative QCD mechanism which can not be 
calculated by using pQCD. 


Note that the potential model calculations 


30| use phenomenological potentials like 


Coulomb potential or Coulomg plus linear potential etc. between color singlet cc by solving 
non-relativistic wave equation to determine the radial wave function at the origin i?(0) for 
the j/p) formation from color singlet cc pair. However, strictly speaking, unlike QED, the 
exact form of the potential energy between cc pair is not known in QCD because we do not 
know the exact form of the classical Yang-Mills potential A^°‘{x) yet, although Yang-Mills 


theory was discoverec 
mental data in QCD 


almost 60 years ago. Hence, unlike QED, we depend on the experi- 


30 


3l| to determine the non-perturbative matrix element of the heavy 


quarkonium production from heavy quark and antiquark pair. The non-perturbative matrix 
element for j/p) formation from color singlet cc pair needs to be universal, i. e., it should 
not change from one experiment to other. Hence it is necessary to use the correct dehnition 
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of the non-perturbative matrix element oi j/'ip production from color singlet and spin triplet 
cc pair. If one does not use the correct dehnition of the non-perturbative matrix element of 
j/pj production from color singlet and spin triplet cc pair then one will predict inhnite cross 
section for j/p) production. This can be seen as follows. 

If the factorization theorem holds, then the production cross section for j/ip at transverse 
momentum Pt can be written as the product of perturbative cross section of cc production 
times the universal non-perturbative matrix element, 

dCA+B^j/i^+XiPT) = do-A+B^ccl^Sij+XiPT) < 0|C>j/p|0 > (6) 

where daA+B^cci^Si]+x{PT) is the perturbative cross section for cc pair production in color sin¬ 
glet and spin triplet state and the non-perturbative matrix element < > represents 

the probability of cc pair in color singlet and spin triplet state to produce j/ip. In eq. ([6]) the 
notation ^Si stands for the spectral notation of the j /%p where S' = 1,L = 0, J=1 

are the spin, orbital angular momentum and total angular momentum respectively. 

The dehnition of the non-perturbative matrix element < 0|Clj/p|0 > can be found as 
follows. The amplitude for the cc to form the specihc hadron H plus other hadrons X is 
given by 

<H + X\cc>. (7) 

The amplitude in eq. o gives the probability of the hadron H production from the cc by 
using the non-perturbative matrix element 

< 0|C1//|0 >= ^ < cc|i7 -|- X >< H + X\cc >= ^ < cc\a\j\X >< Xla^lcc > 
x X 

=< cc\a^fjaH\cc > ( 8 ) 

where is the creation operator of the hadron. Writing the Dirac wave function of the 
quark in terms of creation and annihilation operators a^, a and then using the dehnition 
of the charm quark state |c > formed from the QCD vacuum |0 > by using the creation 
operator ah 

|c >= a^|0 > (9) 

we hnd from eq. ([8]) that the dehnition of the non-perturbative matrix element at the origin 

for the T]c production from color singlet cc pair is given by 

< 0\Oh\ 0 >=< 0\x^^pa^H(^HpJ^x\0 > (10) 
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where x are two component Dirac spinors of the heavy quark wave function. Similar to 
the derivation of eq. flTOj) we hnd that the dehnition of the non-perturbative matrix element 
at the origin for the j/'ip production from cc pair in color singlet and spin triplet state is 
given by 

< 0|C>h|0 >=< > (11) 

where cr* is the Pauli spin matrix. 

The wave function at the origin i?(0) for the j/pj production from cc pair in color singlet 
mechanism is given by 

1 ^( 0 )^ = > • ( 12 ) 

Note that in the perturbative cross section daA+B^cc[^Si]+x{PT) ^q. ([6]) the infrared 
divergences can occur due to the interaction of the c and/or c with nearby light quark or 
gluon. This infrared divergences occur due to the soft (infrared) gluon interaction between 
the c and/or c with nearby light quark or gluon. Hence it is necessary to prove that any 
non-canceling infrared divergences in the perturbative cross section daA+B^cc[^Si]+x{PT) 
eq. dHD cancels with the corresponding infrared divergences in the non-perturbative matrix 
element < 0|(Pj/p|0 > in eq. ([6]). If such a cancelation does not happen then one will predict 
inhnite cross section daA+B^j/'>p+x(PT) j/t/ production from eq. ([6]). Such cancelation 
of infrared divergences is called factorization. Hence it is necessary to prove factorization 
theorem oi j/p) production at high energy colliders. 

Similar to eq. ([9]) in vacuum, using the dehnition of the |c > state in non-equilibrium by 
using the creation operator ah 

|c >= a^\in > (13) 

we hnd from eq. (IH]) that the non-perturbative matrix element for rjc production from color 
singlet cc pair in non-equilibrium QCD is given by 

< in\OH\in >=< > (14) 

where |m > is the initial state of the non-equilibrium QCD medium. Eq. (flTD is similar to 
eq. flTOl) except that the vacuum expectation is replaced by medium average. Similar to the 
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derivation of eq. ffT^ we find that the non-perturbative matrix element for j/'ip production 
from cc pair in color singlet and spin triplet state in non-equilibrium QCD is given by 

< in\OH\in >=< in\x^a'^'tpa}^aH'<P ^> (15) 

which is similar to eq. flTT]) except that the vacuum expectation is replaced by medium 
average. 

The infrared divergences can occur due to the (soft) gluon exchange between the c and/or 
c with nearby light quark or gluon in non-equilibrium QCD medium at RHIC and LHC. 
Hence it is necessary to prove that any non-canceling infrared divergences due to these soft 
gluons exchange in non-equilibrium QCD medium cancel with the corresponding infrared di¬ 
vergences in the non-perturbative matrix element < in\Oji^\in >. If such a cancelation does 
not happen then one will predict inhnite cross section iox jjip production at RHIC and LHC. 
Such cancelation of infrared divergences is called factorization in non-equilibrium QCD. 
Hence it is necessary to prove factorization theorem of j/ip production in non-equilibrium 
QCD at RHIC and LHC. 

In this paper we will prove that the infrared divergences due to the soft gluons exchange 
between charm quark and the nearby light-like quark (or gluon) in non-equilibrium QCD 
at RHIC and LHC exactly cancel with the corresponding infrared divergences due to the 
soft gluons exchange between anticharm quark and the same nearby light-like quark (or 
gluon) at all order in coupling constant in the j/ip production from color singlet cc pair 
in non-equilibrium QCD. This proves the factorization theorem of j /ip production in non- 
equilibrium QCD at RHIC and LHC at all order in coupling constant. 


IV. INFRARED DIVERGENCE AND UNPHYSICAL LONGITUDINAL POLAR¬ 
IZATION OF THE GAUGE FIELD IN QUANTUM FIELD THEORY 


Proof of factorization theorem is usually given in diagrammatic method by using pQCD 


20 


24| . However, since the matrix element < j/t/ production 


from cc pair in eq. flTTll is the non-perturbative matrix element it is useful to use the path 


integral formulation of QCD to prove factorization theorem. In fact, t 
tion theorem using path integral formulation is enormously simplihed 


le proof of factoriza- 


28| in comparison to 


diagrammatic method. Also path integral method naturally proves factorization theorem 
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at all order in coupling constant whereas the diagrammatic method using pQCD suffers 


difficulties beyond certain order of coupling constant calculation 


2 ^. 


Below we will briefly describe the general technique used in the path integral formulation 
in any quantum field theory to prove factorization theorem. Since we are interested in the 
soft-gluons exchange between the charm (and/or anticharm) quark with the nearby light-like 
quark and/or gluon we need to study the infrared behavior of the non-perturbative correla¬ 
tion function of the type < 0|^(a;i)d'(a;i)^(a;2)'h(a;2)|0 > in QCD in the presence of light-like 


Wilson line, in order to study j/'ll) production from color single 


charge produces pure gauge field in classical mechanics 


20 


32 


cc pair. Since a light-like 


33| and in quantum field 


theory (see the next section), one can study the infrared behavior of the non-perturbative 
correlation function of the type < 0|^(a;i)d'(a;i)^(a;2)tl/(a;2)|0 > in QCD by using the path 
integral formulation of the background field method of QCD in the presence of pure gauge 
background field in order to prove factorization infrared divergences at all order in coupling 
constant 


Note that the pure gauge field in quantum field theory corresponds to longitudinal po¬ 
larization of the gauge field. In this section we will show how the longitudinal polarization 
of the gauge field is used in quantum field theory to describe soft (infrared) divergences. 
Since the argument is valid in any quantum field theory let us see how the pure gauge field 
is used in QED to prove factorization of infrared divergences in QED. Consider an incoming 


electron of four momentum and mass m emitting a real photon of 


bur momentum in 


QED. The corresponding Feynman diagram contribution is given by 341 

1 


M = 


where we write 


— 'yuk'^ — 'rn 




p ■ k 


2p ■ k 


(16) 


M 


eikonal — 


p ■ e{k) 
p ■ k 


u(p) 


(17) 


and 


-Mnon-eikonal = -^- m(p)- 


(18) 


From eq. (3.2) of [3^ we write the gauge field as 

,nk) = le'‘(k) - = y„^Jk) + 4„Jk) 


p ■ k 


p ■ k 


(19) 
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where 


( 20 ) 

is the physical gauge held [corresponding to transverse polarization of the gauge held] and 

( 21 ) 

is the pure gauge held [corresponding to longitudinal polarization of the gauge held]. 

Now using eq. flT^ in flThl) we hnd that the total contribution of the Feynman diagram 
is given by 

■ 2 ^ -^^eikonal “ 1 “ -^^non—eikonal (22) 


where 




eikonal 


P ' ^phys (^) / \ P ^pure {k) 


and 


M non—pi 


non—eikonal 


p ■ k 
k'''^^u7,^(^phys{k) 


■u{p) - 


p ■ k 


u{p) = 


P ■ *ipure(^) 
p ■ k 


«(p) 


(23) 


.,(p) + (24) 


2p ■ k 2p ■ k 2p ■ k 

Hence in the soft photon limit (fco, ^ 2 , ^ 3 ) —t 0 we hnd from the eqs. (1T6|1 and (1^ 
that 


-M 


eikonal 


p ■ e{k) 
p ■ k 


u{p) 


P ■ Cpure(^) 
p ■ k 


u{p) —>■ CX) 


as (fco, ^ 1 , ^ 2 , ^ 3 )-t 0(25) 


which implies that the physical gauge held [corresponding to transverse polarization] does 
not contribute to the soft (infrared) divergences in quantum held theory and the soft (in¬ 
frared) divergences can be calculated by using pure gauge held [corresponding to longitudinal 
polarization] in quantum held theory. 

Similarly from eqs. flTB]) and we hnd that 


M non—I 


non—eikonal 


_ k^lul^,e^{.k) _ fc^7i27/4ephys(fe) 


2p ■ k 


-u{p) = 


2p ■ k 


u{p) hnite 


as (fco, A:i,/c 2 , ^ 3 )-)■ 0 

(26) 


which contribute to the hnite part of the cross section which implies that pure gauge held 
[corresponding to longitudinal polarization] does not contribute to the hnite cross section 
and the hnite cross section can be calculated by using physical gauge held [corresponding 
to transverse polarization]. 
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Hence we find that the non-eikonal-line part of the diagram as given by eq. (1261) is 
necessary if we are calculating the hnite value of the cross section but it is not necessary if 
we are calculating the relevant infrared divergence behavior. The relevant infrared divergence 
behavior can be calculated by using the eikonal-line part of the diagram as given by eq. (12^ . 

Hence we hnd that we do not need to calculate the hnite value of the cross section (or the 
full cross section) [which will require the non-eikonal-line part of the diagram as given by eq. 
fl26|l ] to study the relevant infrared divergence behavior. The relevant infrared divergence 
behavior can be calculated by using eikonal approximation as given by eq. fl2^ . 

From eq. fl2B]) we hnd that 

Aii:r-£s.“ = = 0. (27) 

We are interested in the infrared divergence behavior due to the presence of the light-like 
Wilson line. We will show in the next section that the eikonal current of the light-like charge 
generates pure gauge held in quantum held theory. Hence from eqs. (l25jl and (ITTIl we hnd 
that the soft (infrared) divergence behavior due to the presence light-like Wilson line can 
be studied by using pure gauge held in quantum held theory without modifying the hnite 
value of the cross section 


V. GAUGE FIELD GENERATED BY EIKONAL CURRENT OF LIGHT-LIKE 
WILSON LINE IN QUANTUM FIELD THEORY 


In order to study factorization of infrared divergences by using the background held 
method of QED, the soft photon cloud traversed by the electron is represented by the pure 


gauge background held 


25| due to the presence of the light-like Wilson line. As mentioned 


above, in classical mechanics the assertion that the gauge held that is produced by a highly 
relativistic (light-like) charge is a pure gauge held at all time-space position except at the 


position transverse to the motion of the charge at the time of closest approach 20, l32|, l33 |. 


One may ask a question if this assertion is correct in quantum held theory. In this section 
we will show that this assertion is correct in quantum held theory. We will use path integral 
formulation of the quantum held theory to show this. 

The generating functional for the gauge held in the quantum held theory in the presence 
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of external source J^{x) in the path integral formulation is given by 

Z[J] = J (28) 

where 




The effective action Seff[J] is given by [351] 


= F>^^[A]F,M]- 


< 0|0 >J= 


m 


where 


Seff[J] = —- J d'^xdf'x' J^{x)D^y{x — x')J^{x') 

where D^y{x — x') is the photon propagator. 

The photon propagator in the coordinate space is given by 

D^y{x - x') = d^dy\5^'^\x - x'). 

Using eq. (15^ in (ITT]) we hnd 

Se,AA = -\j + ''-^^d.dArix). 

From the continuity equation we have 


(29) 


(30) 


(31) 


(32) 


(33) 


a^JAx) = 0. (34) 

Using eq. (IMj) in (133|) we hnd 

Seff[J] = -\j d'^xJ>^{x)^Jf,{x). (35) 

First of all, by using the path integral formulation of the quantum held theory we will 
derive Coulomb’s law for static charge. Note that the derivation of the Coulomb’s law 
by using path integral formulation of the quantum held theory is not necessary to prove 
factorization theorem. We have included it here only to demonstrate the correctness of the 
prediction of the path integral formulation in quantum held theory which we will use (see 
below) to show that the eikonal current of the light-like charge generates pure gauge held in 
quantum held theory. 
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In order to derive Coulomb’s law by using path integral formulation of the quantum held 
theory let us consider the static charge at the position X. The current density for this static 
charge is given by 


^(x) = - X). 


(36) 


m 

,2 


we hnd 


Using eq. 

SeiAA = jj = 


.^ljd*x[X-X 

^ \x — 


X\ 


\x 


X\ 


which gives the effective lagrangian density 

.2 ^ 


^e//((r) = -y[V- 

^ \x — X\ 


• V- 


E^{x) 


(37) 


(38) 


which reproduces the Coulomb’s law. Hence we have shown that the assertion that a charge 
at rest generates a Coulomb gauge held is correct in quantum held theory. 

Similarly using the above procedure in quantum held theory we will show that the asser¬ 
tion that the eikonal current of the light-like charge generates pure gauge held is correct in 
quantum held theory. This can be shown as follows. 

In QED the infrared (or soft) divergence arises only from the emission of a photon 
for which all components of the four-momentum are small. As described in the previous 
section, the Eikonal propagator times the Eikonal vertex for a soft photon with momen¬ 
tum k interacth^ with a light-like electron moving with four momentum is given by 


20 , y, y, y, y, 1361421 






- , . = e —-- (39) 

p ■ k -\- le I ■ k + le 

where E is the four-velocity of the light-like electron. Note that when we say the ” light-like 
electron” we mean the electron that is traveling at its highest speed which is arbitrarily close 
to the speed of light (|/| ~ 1) as it can not travel exactly at speed of light (|/j = 1) because 
it has hnite mass even if the mass of the electron is very small. From eq. fl3^ we hnd 
d'^k l-A{k) ■ f 

— ‘ ■ A{k) = le / dXl ■ A{IX) ( 40 ) 


/ N. , , = —ei I dX I -—- 

( 27 r)"‘ I ■ k + ie Jo J ( 27 r)' 


where the photon held A'^{x) and its Fourier transform A'^{k) are related by 

r d‘^k 


Ak-x 


(41) 
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From eq. fHOj) we find 


POO P 

ie J dXl ■ A{IX) = ^ J d‘^xJ^{x)A^{x) 


(42) 


where the eikonal current density J^{x) of the light-like charge e is given by 


J^^[x) = eP J dX 6^*\x-lX), 


f = = 0 


(43) 


Hence by using the path integral formulation of the quantum held theory we hud by using 
eq. fl43|) in fl3^ that the effective action is given by 

.2 . . 1 


Seff[J] = j J dX S^^\x J dX' 6^^\x - IX') 
2 

= J d^x[d^^- 




which gives the effective lagrangian density 


d^effix) — I 




From eq. (H5ll we hud the effective lagrangian density 

e2 (/2)2 


d^effix) — 


2 (/•a;)4’ 


for 


/ • X 7^ 0 


(44) 


(45) 


(46) 


which gives 


d^effix) = 0 (47) 

at all time-space position except at the position transverse to the motion of the charge 
(/ • X = 0) at the time of closest approach (xq = 0) where we have used 

f = = 0 (48) 

where ^ with = (c, Fc) where and c is the speed of light. In this paper we 

have used natural unit where c = 1. 

Hence from eqs. fl30|) . ([35]), (1471) . (1431) and (145|l we hnd that the eikonal current for 
light-like charge generates pure gauge held in quantum held theory. 
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VI. PROOF OF FACTORIZATION OF J/^' PRODUCTION AT HIGH ENERGY 
COLLIDERS 


at 


in 


Arguments about the proof of factorization oi j/'ip production in color singlet mechanism 
ligh energy colliders by using pQCD diagrammatic method in QCD in vacuum is given 
30| . In this section we will prove the factorization theorem oi j/ip production from color 


singlet and spin triplet cc pair at high energy colliders at all order in coupling constant 
by using path integral formulation of QCD in vacuum. We will extend this path integral 
technique to prove factorization theorem of j/p) production from color singlet and spin 
triplet cc pair in non-equilibrium QCD at RHIC and LHC at all order in coupling constant 
in section VII. 

The generating functional in QCD including the heavy quark is given by 


43 


44| 


Z[j, r]u, r]u, rjd, Vd, Vs, Vs, Vh, Vh] = 


J [dpJi] [dpj2\ [d'ip2] [d^p3\ [dips] [dT] [dT] det( 


6oj^ 




(49) 


where is the quantum gluon held, the symbols Z = 1, 2, 3 = u, d, s stand for three light 
quarks u, d, s and the symbol h stands for heavy quark and 


F‘JQ] = dfQi(x) - dMx) + 9 !‘^Qi(x)Ql(x), 


pa2 
fiiy L 


[Ql = KAQ\F‘“'‘[Q\- (50) 


Note that the determinant det(^^%^^) in eq. (I49|) can be expressed in terms of path 
integration over the ghost helds 4^. However, we will directly work with the determinant 
det(^^%^) in eq. (09]). 

For the heavy quark Dirac held T(x), the non-perturbative matrix element of the type 
< 0|T(a;i)\k(a;i)T(x2)T(a;2)|0 > in QCD is given by 25| 


< 0|4>(a:i)\k(xi)4>(a:2)T(a:2)|0 > 

= J [dQ][d^pl][d^pl][d^p2][d^p2][d^p3][d^p3][d^][d^] T(a;i)T(a;i)T(a;2)T(a;2)det(^^^|^-^) 


where the suppression of the normalization factor Z[0] is understood as it will cancel in the 
hnal result (see eq. ([73])). 
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As mentioned in sections IV and V the infrared divergences can be studied by using 
eikonal approximation in quantum field theory. We are interested in the infrared divergence 
behavior of the non-perturbative heavy quark-antiquark correlation function of the type 
< 0|^(xi)\k(a;i)^(a;2)d'(a;2)|0 > in the presence of light-like Wilson line in QCD. 

The Eikonal propagator times the Eikonal vertex for a gluon with four momentum 


interacting with a light-like quark moving with four velocity is given by 20, l24l-l27l. 


36 


42| 


gr 


(52) 


I ■ k + ie 

Hence for a light-like quark attached to infinite number of gluons we find the eikonal factor 

d^ki d^k2 TV ■ A^{ki)TH ■ A\k2) 


d^k l-A^{k) 2 

(27r)^ I ■ k + ie^ ^ . 


(27r)'^ (271)^ {I ■ ki + ie){l ■ {ki + ^ 2 ) + ie] 


+ 


l + gT<^ j 

POO POO POO 

= 1 + / d\l ■ A^{1\) + g‘^e / d\i / dA2TV ■ A\IXi)TH • A'’(/A2) + ... 

JO JO JAi 

POO cP'i^ 

= 1 + / dXl ■ A’^ilX) + / dXi / dA2TV ■ A^{IXi)TH ■ A\IX2) + ... 

JO 2! Jo Jo 

POO 

= V exph^ / dXl ■ A“(/A)T“] (53) 

Jo 

which describes the infrared divergences arising from the inhnite number of soft gluons 
exchange with the light-like quark where V is the path ordering and the gluon field 
and its Fourier transform A^°‘{k) are related by 

d^k 


A^^'^ix) = 


(2w) 


fV“(J;)e 


ik-x 


(54) 


. As mentioned in sections IV and V the light-like quark traveling with light-like four-velocity 
produces SU(3) pure gauge held at all the time-space position except at the position 
X perpendicular to the direction of motion of the quark (/ • x = 0) at the time of closest 


approach 


20 


32 


33|. When = A^“(AZ) as in eq. fl53|l we hnd l-x = XI ■l = X^0 

which implies that the light-like quark Ends the gluon held in eq. (I5^ as the SU(3) 

pure gauge. The SU(3) pure gauge is given by 


T-A-lx) = -[d^U(x)] U-'(x), 


which gives 


^9 


U{xf) = 


U{x) = 


(55) 


U{xi) = 


(56) 
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Hence when A^°‘{x) = as in eq. fl5^ we find from eq. 

line in QCD for infrared divergences is given by 


/g°°<iAZ-A“(x+iA)T“ ^ ^igT°'ui°'{x) 


that the light-like Wilson 


(57) 


Hence the effect of infrared gluons interaction between the partons and the light-like Wilson 
line in QCD can be studied by putting the partons in the SU(3) pure gauge background 
field which implies that the infrared behavior of the non-perturbative correlation function 
of the type < 'l'(a;)\k(a:')'h(a:")\l'(x'")... > in QCD due to the presence of light-like Wilson 
line in QCD can be studied by using the path integral method of the QCD in the presence 
of SU(3) pure gauge background field as given by eq. fISSD . 


Background field method of QCD was o rig inally formulated by ’t Hooft 45| and later 


extended by Klueberg-Stern and Zuber 


46, 


47l | and by Abbott This is an elegant for¬ 


malism which can be useful to construct gauge invariant non-perturbative green’s functions 
in QCD. This formalism is also useful to study quark and gluon production from classical 



chromo field via Schwinger mechanism to compute (3 function in QCD to perform 


calculations in lattice gauge theories and to study evolution of QCD coupling constant 
in the presence of chromofield 501. 

The generating functional in the background field method of QCD is given by 


Z[A, J, r]u, Vu, Vd, Vd, Vs, Vs, Vh, Vh] = 


i] [#i] [(iV’ 2 ] [dip 2 ] [dV’ 3 ] [dips] [d^ [dd'] det(— 


600 ’^ 


expp 
3 


d^x\--F^l,[A + Q]- + J-Q 


4 - g..-- . -V. 2 a 

-mi+ gT^Y{A + + Viipi + ipiVi 


1=1 


-M + grYiA + + ito'S + iri„]] 


(58) 


where the gauge fixing term is given by 


G“(Q) = 

which depends on the background field A^°‘{x) and 

KM + Ql = d,[Al + Qa - dMl + «1 + + Qiu: + Ql], 


(59) 


(60) 


We have followed the notations of 


m 


m. 


46l | and accordingly we have denoted the quantum 


gluon field by Q^°‘ and the background field by The gauge fixing term in 
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( 61 ) 


eq. fl58|) [where G°'{Q) is given by eq. fl5^ ] is invariant for gange transformation of ^4“: 


5/1“ = (type I transformation) 


provided one also performs a homogeneous transformation of Q“ 


43 


46|: 


^q; = gf^Qy- 


(62) 


The gauge transformation of background held ^4“ as given by eq. fl6T|) along with the 
homogeneous transformation of Q“ in eq. fl6^ gives 

5(kl“ + g“) = + d.cj- (63) 


which leaves —^F°‘‘^^j^[A + Q] invariant in eq. 


For the heavy quark Dirac held T(a;), the non-perturbative matrix element of th e ty pe 
< 0|'F(xi)\k(a;i)'F(a;2)4/(x2)|0 > in the background held method of QCD is given by 25| 


< 0|T(xi)T(xi)T(x2)T(x 2)|0 >2l 

= J [dg][dV)i][#i][#2][#2][#3][#3]M^][c^4^] ^(a;i)T(xi)T(x2)T(x2) det( ^^J^^ ) 

exp[s I d*xl~FyiA + <31 - ^{G‘(Q)r + + gT‘Y{A + Q)y, 

+*[i7'‘ap - M + 9T“7'‘(A + Q);]®]] (64) 

where the suppression of the normalization factor Z[0] is understood as it will cancel in the 
hnal result (see eq. fTfdll ). By changing Q y- Q — A in eq. fl64ll we hnd that 


< 0|T(xi)T(xi)T(x2)T(x2)|0 >a 

= J [dg][#i][#i][#2][#2][#3][#3]M4'][ciT] T(xi)T(xi)T(x 2 )T(x 2 ) det( ) 

exp[i J d‘^x[~F‘^l^[Q] - ^iG}{Q)f + Y. - mi + c/TVg“]^z 

-M + (7T“7'^g“]T]] (65) 

where the gauge hxing term from eq. fl5^ becomes 

G%Q) = + gr^'^AlQ^^^ - d^A^^ = D^[A]Q^^ - d^A^^ (66) 

and from eq. fl6^ [by using eq. fl6T]) ] we hnd 

SQ‘ = (67) 
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For finite transformation eq. fl67|) becomes 

r“Q;“(x) = f/(a;)T“Q“(a:)f/-'(x) + -[d^U{x)]U-\x), U{x) = (68) 

'^9 

The fermion fields transform as 




^'(x) = 


(69) 


Changing the variables of integration from unprimed to primed variables in eq. fl65l) we hnd 


< 0|T(a;i)T(a;i)T(a;2)T(a;2)|0 >a 




'] Wi] Wi] W 2 ] W 2 ] ^'(a:i)d''(a;i)T'(a; 2 )T'(a; 2 ) det( ' ) 


1 


exp[i J - ^(G/(«'))" + S + aT°YQ'‘W 

+'¥\iYa^ - M + gT‘YQ“W]] 


(70) 


because the change of variables from unprimed to primed variables does not change the 
value of the integration. 

From eqs. (|6^ and (|6^ we hnd 


[dQ'] = [dQ], = [#l][#l], [#2][#2] = ['^V’2][#2], [#3][#3] = [#3][#3], 

[dT'][dT'] = [dT][dT], -mi + gT'^^Q’^Wi = -mi + 


- M + c/TVQ^T]^' = -M + c/T“7^Q“]T, 


flU 


[Ql = KM 


{g%q')Y = (a,Q'^“(x))^ 

Using eq. fl7T]l in flTOjl we hnd 

< 0|T(a;i)T(a;i)T(a;2)d'(a;2)|0 > 




6u^ 


6u^ 


(71) 


J [dtlai] [di/’2] [d'ip2] [di/’s] MV's] [dd^] [dT] T(a;i)T(a;i)T(a;2)T(a;2)det(^^%^-^) 




7 f d-ix[-jF-f,JQ]-£(d^Qn^+Y;^^^l^l[i'yao^-mi+gT-rQ-]^lJ+^[i'y>^d^-M+gT-'y<^Q-m 


From eqs. (1721) and (I5T1) we hnd 


(72) 


< 0 |^(Xi)'If(xi)W(x 2 )^(x 2)|0 >^< 0 |^(xi)'If(xi)^(x 2 )'If(x 2)|0 >A 


(73) 


which proves factorization of infrared divergences at all order in coupling constant. Eq. 
fl73|) is valid in covariant gauge, in light-cone gauge, in general axi al g auges, in general 


non-covariant gauges and in general Coulomb gauge etc. respectively 


28|. 
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From eq. fl73|) we find that the non-perturbative matrix element for j/'ll) production from 
cc pair in color singlet and spin triplet state in eq. (JO as given by 

< 0|C>h|0 >=< > (74) 

is consistent with factorization theorem of j/'ip production at high energy colliders at all 
order in coupling constant. 

Eq. (|73il proves the factorization theorem of j/'ip production at high energy col¬ 
liders at all order in coupling constant where < 0|'F(a:i)\l/(a;i)^(a:2)4/(a^2)|0 > is the 
heavy quark-antiquark gauge invariant non-perturbative correlation function in QCD and 
< 0 |^(a;i)'k(a;i)^(a; 2 )tl/(a; 2)|0 >a is the corresponding heavy quark-antiquark gauge invari¬ 
ant non-perturbative correlation function in QCD in the presence of light-like Wilson line. 
From eq. fl75D one hnds that the infrared divergences due to the soft gluons exchange be¬ 
tween charm quark and the nearby light-like quark (or gluon) at high energy colliders exactly 
cancel with the corresponding infrared divergences due to the soft gluons exchange between 
anticharm quark and the same nearby light-like quark (or gluon) at all order in coupling 
constant in the j/'ip production from color singlet cc pair at high energy colliders. This 
proves the factorization theorem of j/'ip production at high energy colliders at all order in 
coupling constant. 

VII. PROOF OF FACTORIZATION OF J/^ PRODUCTION IN NON¬ 
EQUILIBRIUM QCD AT RHIC AND LHC 

The generating functional in non-equilibrium QCD (including heavy quark) in the path 
integral formulation is given by eq. (jS]). The nonequilibrium-nonperturbative heavy quark- 
antiquark correlation function of the type < m|Tr(a;i)Tr(a:i)Ts(a; 2 )Ts(x 2 )|m > in QCD is 
given by eq. ([5]) which gives 

< m|'Fr(a;i)Tr(xi)Ts(a;2)4/s(a;2)|m >= 

= J [dQ+] [dQ_] [#i+] [d'ipi_] [#i+] [d'ipi_] [d'ip 2 +] [# 2 -] [dip 2 +] [d'ip 2 -] [dip3+] [dips-] [d'ip3+] [dips-] 

df) 0 ^^“ 

[dT+][dT_][dT+][dT_] Tr(a:i)Tr(a;i)T^(a: 2 )Ts(a: 2 ) x det( ) x det( ) 

exp[* / d'xl-pF-llQ^] - f^Uq-]) - 
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+ -mi + - II -m + gT^j^Q^^_]'ipi_ 

1=1 1=1 

+^+[trd^ -M + gT-YQU'^+ - -M + gT^^YQU^-]] 

X < Q+,r+,^i,K^^lr+,^+,^+,^+,o\ P io,^_,'i'_,c,c,^i,V'i,c,c,Q- > 


(75) 


where the repeated (closed-time path) indices r, s = -f, — are not summed and the suppres¬ 
sion of the normalization factor Z[0] is understood as it will cancel in the hnal result, see 
eq. ([82]). 

The generating functional in the background held method of QCD in non-equilibrium 
QCD (including heavy quark) in the path integral formulation is given by 45-T7,|^, 54| 


Z[p, A, J+, J_, , 77+, gl, r/t 7-, h+, h-, h+, 7-, h+, h+, V-] 




= j [dQ+] [dQJ\ [dil)i+] [#i-] [dipi+] [dipiJ] [#2+] [di^ 2 -] [di^ 2 +] [#2-] [^7/^3+] [#3+] [c?^3-] 

[(iT+][dT_][dT+][dT_] x det( '^^^^^+^ ) x det( ^^^^^~^ ) x exp[i J d'^x[-^{F'^l^[Q+ + A+] 

-F^lAQ- + 7l_]) - - {G\Q.)f) + - m, + gT<^Y{Q+ + 7l+)“]77z+ 

/a 

- -mi + gT^YiQ- + + ^+[^^9, -M + ^TV(Q+ + 

1=1 

-^-[^^9, -M + gT^riQ- + + J+ ■ Q+ - J_ • g_ 

3 

+ Y^[f]i+ ■ ipi+ - f]i_ ■ iIji^ + ipi+ ■ gi+ - iIji^ ■ gi_] + %+ ■ T+ - f]h- ■ -h T+ ■ gh+ - ■ Vh-]] 

1=1 




(76) 


where the gauge hxing term G°'{Q) is given by eq. (|59|) which depends on the background 
held y4^“((r), the F“j,[g -l- A] is given by eq. dHOj), the 5^“ and are given by eqs. fl^ 
and fIM]) respectively. 

The nonequilibrium-nonperturbative heavy quark-antiquark correlation function of the 
type < m|Tr(a;i)\['r(a;i)Ts( 0 : 2 )ft's( 0 : 2 ) |m > in the background held method of QCD is given 
by 


< m|Tr(xi)Tr(xi)Ts(x2)T5(x2)|m >a 

= j[dQ+] [dQ-] [# 1 +] [di/ji-] [dT/'i+j [dipi-] [# 2 +] [# 2 -] [d'ip 2 +] [d7/'2-] [dip3+] [# 3 -] [d'ip3+] [di>3-] 
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[d^+] [d^.] [d^+] [d^.] 

X X exp[i jd‘x\~(F‘‘lJfi+ + A+] 

-F^llQ. + /!_]) - ^((G»(QJ)2 - (G“(Q_))2) + j:i;+liYd, - m, + gr^Q^ + 

/a 

-mi + gT^riQ- + ^-)“]^^- + ^+[^^9, -M + gT^riQ+ + ^+)“]^+ 

1=1 

-M + ^TV(Q- + x4-)“]^-]] 

X < g++A+, r+, <,^+> ^+> oi p io, ^-, c,c,c, q-+ x4_ > 

(77) 

where the suppression of the normalization factor Z[0] is understood as it will cancel in the 
hnal result, see eq. fl8^ . 

In order to study the infrared behavior of the nonequilibrium-nonperturbative heavy 
quark-antiquark correlation function we proceed as follows. By changing the integration 
variable Q ^ Q — A in the right hand side of eq. fl77|) we hnd 

< in\^rixi)'^rixi)^six2)'^six2)\in >A 

= j[dQ+] [dQ_] [d^i+] [d^i-] [#i+] [dipi-] [d^/’2+] [# 2 -] [# 2 +] [d^ 2 -] [d^3+] [dips-] [dip3+] [dips-] 

[d^+] [d^_][d^+][d^_] ^,(xi)d^,(a:i)^,(x2)d^s(a:2) 
xdet(^^^^) X det(^S^^) x exp[i J d*x[~{F^l^[Q+] 

-F“j,[Q_]) - ^((g^(q+))2 - {G%Q.)f) + j2P’i4zrd, -mi+ 

/a 

-j2iPi_[tYd, -mi + gT‘^rQl_]iPi. + ^+[ 1 ^ 8 , -M + 

1=1 

-^_[^ 7 '^a^ - M + ^ 7 TVQ“-]^-]] 

X p \o,^.,^.,r_,r-,'9ti^tip-,ip-,Q- > 

(78) 

where GJ{Q) is given by eq. fl66|) and eq. fl6^ [by using eq. (16T]) ] becomes eq. fl67)) which 
for hnite transformation gives eq. flB5]l . Changing the integration variable from unprimed 
variable to primed variable we hnd from eq. fl78p 

< m|'hr(xi)d'r(a:i)f s(a;2)^s(a:2)|w >a 

= J [dg+] [dQ'_] [dip[^] [dip[_] [dip[^] [dip[_] [#^+] [dip'^.] [dip'^^] [# 2 -] [# 3 +] [d'9s+] Ws-] 
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xdet( 


6G%Q' 


X det("^^ 


6G%Q') 


) X exp[z J d^x[-\{F^l[Q'^] - F^l[Q'_]) - ^{{G%Q'^)f 


5uj^ 6ut 

3 3 

-(G‘;(Q'_)f) + J2i'',4ird, - ra, + gr','‘Q'‘M+ “ ECln'S, - m, + 

i=l 1=1 

-M + gT-YQ'^^]K - -M + gT‘^YQ'^_]^-]] 

X < g'+, iP':, iP’:, iP'^, iP'l iPt, n, n, 0| p |0, v[/'_, ^'_,ip'P, ip'P, Pp't , V'-, Q'- > • 


(79) 


This is because a change of integration variable from unprimed variable to primed 
variable does not change the value of the integration. Note that since we 
are working in the frozen ghost formalism at the initial time 


53 


54j | the < 


Q+, ipl, ipl, ip^, ip^, ipp, ipp, T+, T+, 0| p |0, T_, T_, ipi, ipi, ipi, ipi, V’-, , Q- > in eq. ([3]) 

corresponding to initial density of state in non-equilibrium QCD is gauge invariant by deh- 
nition. Hence from eqs. flbSp and fl69p we hnd 


< O'+W-r. c. '/>?, i’t i’t i’t «'+.»'+. 0I p I0. V- ■ Ct i’-, V-, Q'- > 
=<«+,«,■«W't'AtV'l,4'+,0| p |0,> 


( 80 ) 


We follow the similar procedure that was employed to prove factorization theorem of j/ip 
production in QCD in vacuum in the previous section. When background held A^°'{x) is the 
SU(3) pure gauge as given by eq. (jSS]) then by using eqs. (ITip and (IHUP in eq. (ITUP we hnd 

< in\^r{xi)'^rixi)^six2)'^s{x2)\in >A= 

= J[dQ+] [dQ_] [# 1 +] [# 1 -] [dipi+l [dipi^] [# 2 +] [dip 2 -] [d'ip 2 +] [# 2 -] [# 3 +] [d'lps-] [# 3 +] [# 3 -] 
[dT+][dT_][dT+][dT_] Tr(xi)T^(xi)T^(x 2 )T,(a; 2 ) x det( ) x det( ) 

exp[* j d'x[-\{F'‘%[Q^.] - F“=„[Q_]) - T((5„Q7)= - 

3 3 

+ ^-mi + gT‘^YQU]'^i+ ~ ^ P’i-[iYd^ - rrn + pTVg))-]V^z- 

1=1 1=1 

+^4t'y^^d^ -M + pTVg;)+]^+ - -M + pTVQ;)-]^-]] 

( 81 ) 
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From eqs. fl75|) and flHT]) we find 


< in\'^r{xi)'^r{xi)^six2)'^six2)\in >=< m|\l>r(a;i)\l>r(a;i)\l>s(a;2)\l/s(a;2|m >a (82) 


which proves factorization of infrared divergences at all order in conpling constant in non- 
eqnilibrium QCD. Eq. flH^ is valid in covariant gange, in light-cone gange, in general axial 
ganges, in general non-covariant ganges and in general Conlomb gange etc. respectively |28|. 

From eq. (|82|l we find that the noneqnilibrinm-nonperturbative matrix element for j/'ip 
prodnction from cc pair in color singlet and spin triplet state in non-eqnilibrinm QCD in eq. 
flT^ as given by 


< in\OH\in >=< > 




(83) 


is consistent with factorization theorem oi j/pj prodnction in non-eqnilibrinm QCD at RHIC 
and LHC at all order in conpling constant. 

Eq. (I8^ proves the factorization theorem oi j/pj prodnction in non-eqnilibrinm QCD at 
RHIC and LHC at all order in conpling constant where < m|'Fr((Z(i)d'r((Ci)'Fs((Z( 2 )d' 8 ;((Z( 2 )|w > 
is the heavy qnark-antiqnark gange invariant noneqnilibrium-nonpertnrbative correlation 
fnnction in QCD and < >a is the corresponding heavy 

qnark-antiqnark gange invariant noneqnilibrinm-nonpertnrbative correlation fnnction in 
QCD in the presence of light-like Wilson line. From eq. we find that the infrared 

divergences dne to the soft glnons exchange between charm quark and the nearby light-like 
quark (or gluon) in non-equilibrium QCD at RHIC and LHC exactly cancel with the cor¬ 
responding infrared divergences due to the soft gluons exchange between anticharm quark 
and the same nearby light-like quark (or gluon) at all order in coupling constant in the j/p) 
production from color singlet cc pair in non-equilibrium QCD. This proves the factoriza¬ 
tion theorem oi j/p) production in non-equilibrium QCD at RHIC and LHC at all order in 
coupling constant. 


VIII. CONCLUSIONS 

J/ip suppression/production is one of the main signature of quark-gluon plasma detection 
at RHIC and LHC. In order to study j/pj suppression/production in high energy heavy-ion 
collisions at RHIC and LHC, one needs to prove the factorization theorem of j/ip production 
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in non-equilibrium QCD medium, otherwise one will predict infinite cross section of j/'ip- 
In this paper we have proved factorization theorem of j/p) production in non-equilibrium 
QCD at RHIC and LHC at all order in coupling constant. This proof is necessary to 


study the detection of quark-gluon plasma 
suppression/production as its signature 


Q,Q, 


5 si. ISG ] at RHIC and LHC by using j/'ip 
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